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Abstract 


In this paper we propose an explicit and positivity preserving scheme for 
the mean reverting CEV model which converges in the mean square sense with 
convergence order a (a — 1/2). 
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1 Introduction 

Let (12, T, P, J r t) be a complete probability space with a filtration and let a Wiener 
process (Wt)t> o defined on this space. We consider here the mean reverting CEV 
process, 



( 1 ) 


where k,l,cr > 0 and a G (1/2,1). It is well known that this sde has a unique 
strong solution which is strictly positive. Our starting point was the paper of [I] in 
which the author proposes an implicit and positivity preserving numerical scheme 
to approximate the above process. This stochastic process plays important role in 
financial mathematics. If one wants to use the above model to price complicate path- 
dependent options maybe it is useful to approximate it numerically. The usual Euler 
scheme (see 0) does not preserve positivity and therefore numerical schemes with 
this feature are needed. In this direction on can see [TJ and [3]. Our goal here is to 
present an explicit, positivity preserving numerical scheme that converges in the mean 


1 


square sense to the true solution with, at least, a(a —1/2) order of convergence. There 
must be extensive numerical experiments to compare all these methods and decide 
which of them is the best in each set of parameters. 

We will construct our scheme using the semi discrete method that we have pro¬ 
posed in [5] (and further extended in [6j, [7]). Let 0 = t 0 < ti < ... < t n = T and set 
A — A Consider the following stochastic process 


yt = 


n(l - a)(W t - W tk ) + y tk ( 1 - kA) + A (kl 




1 —a 


— r, 1 -a 


= \z t 


( 2 ) 


for t G (t k ,t k+ i\ where 


= ^(1 - a)(W t - W tk ) + y tk ( 1 - kA) + A (kl 


acj2 y 2 tk 1 


1 —a 


for t G (t kl t k+ 1 ]. To show that this stochastic process is well defined we will check 
whether 

aa 2 2a—1 

y tk ( 1 - kA) + A (kl -) > 0. 


For kl > A < —— ~r and noting that 0 < 2a — 1 < 1 we have 
fe +~2- 

cia 2 y 2a ~ 1 \ ( ao 2 v 2a ~ l \ 

yt k ( 1 - kA) + A (kl - f—)j I {ytk>l} + ^ fc (l - kA) + A (kl -)J I {atfc <i} > 

hy tk >i}y 2 r 


1 — kA — 


Aaa z 


+ A l kl- 


acr 


> 0 


Therefore we impose the following assumption. 

Assumption A We assume that Xq G R + . Moreover, we suppose that 

r , a(j2 a 2 

kl >-, A < —- 

2 2 k + acr 2 

This stochastic process, using Ito’s formula, has the differential form, 


yt = Vt k {1 ~ kA) + A (kl 




) + 


rt acr 2 


y 2a 1 ds + cr y a s sgn(z s )dW s , t e (t k ,t k+ 1 ], 




' tk 


Concluding, the numerical scheme that we propose for the mean reverting CEV model 
is the following, 


Vt k +1 


m - aXWfcH - W lh ) + y tt (l - kA) + A (kl - 


9 2n — '\ \ 1 

ao 2 yf k x 
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2 Main results 

Lemma 1 Under Assumption A we have the estimate for the following probability, 


A “"2 

P (z t < 0) < C ——. 

g £±‘2a — 1 

Therefore, the above probability tend to zero faster than any power of A. 
Proof. Indeed, we have 


P (zt < 0) = P [W t - W tk+1 < 


(y tk (l - kA) + Ad) 
er(l — a) 


1—a 


< 


V2tt 


v 

e 2 


CA 1 / 2-0 


dy 


< C 


£a-l/2 


eA 


C 5 
2a —1 


where d = kl — To obtain the last inequality one can use problem 9.22, p.112 of 

0. □ 


Next we will use the following compact form for our scheme, 

U U na 2 

Vt = x o + / (kl - ky s )ds + / — (i/f 1 ' 1 ~ y] a r )ds 
Jo Jo z 

kl - ky tk - ds + a y a s sgn(z s )dW s , t e (t k , t k+1 ]. 

Consider the following process, 

U aa 2 rt 

v t = x 0 + Tkl+ / —y] a i ds+ / ay a s sgn(z s )dW s , 

Jo ^ Jo 

Then it is clear that 0 < yt < Vt- We will show that vt has bounded moments and 
therefore yt has also bounded moments. 

Lemma 2 (Moment bounds) Under Assumption A we have the moment bounds, 

E y 2 + Ex 2 < C, 



for some C > 0 

Proof. Consider the stopping time Or = inf{t > 0 : Vt > R}. Using Ito’s formula on 
v t/\e R we obtain, 


v tAe R — ( x o + Tkl ) 2 + 


(^-v sA 0 R y 2 s “ d ^ + y y 2 s f e ^j ds + 2a v sA e R yf A g R sgn(z sA g R )dW s . 
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Taking expectations on both sides and noting that yt < Vt, we arrive at 

E<A«* < E(x„ + Tklf + C f (Ev 2 MR ) a d3 

Jo 

Using now a Gronwall type theorem (see PS. Theorem 1, p. 360), we arrive at 

^ V tA9 R — C (3) 

But Kv^ 6r = + R p P(6 r < t). That means that P(t A 0 R < t) = 

P( 6 r < t) —> 0 as R —> oo so t A 0 R —> t in probability and noting that d R increases 
as R increases we have that t A d R —> t almost surely too, as R —> oo. Going back to 
(3) and using Fatou’s lemma we obtain, 

En t 2 < C 


The same holds for x t . 

Next we define the process 


h t =x o + / (hi - ky s )ds + a y a s sgn{z s )dW s 


We will show that ht,yt remain close. 


□ 


Lemma 3 We have the following estimates, 


E|y s - ys\ 

E|h s - y s | 2 
E| h s - y s | 2 
E|h s | 2 


< GA for any s G [0, T] 

< GA 2a_1 for any s G [0, T\ 

< GA 2a_1 for any s G [0, T} 

< A, for any s G [0 ,T]. 


Proof. Using the moment bound for y t we easily obtain the fact that 

Ebs - Vs \ 2 < GA 

and then 


E\h s — y s \ 2 < GA 2a_1 . 


Next, we have 

E|h s - y tk f < 2E| h s - y s | 2 + 2E| y s - y tk | 2 < GA 2a_1 . 

Finally, to get the moment bound for h t we just use the fact that is close to y t , 
i.e. 


Eh 2 < 2E| h t — yf 2 + 2Ey 2 < G. 


□ 
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Lemma 4 (Inverse Exponential Moments) For the true solution xt it holds the 
following bound, 


E 



C 


2(l-a) • 


< OO, 


for any C > 0. 

Proof. We first transform our equation with z t = xf' 1 ^ a \ Then using Ito’s formula 
we deduce that 

zt = z 0 + f ( ff 2 (l - a)(l — 2a) - 2(1 - a)kz s + — 2a _f^ J ds + 2(1 - a)a f y/zfdW s , 
Jo \ z 2 ~ 2a J [) 

and if we denote by b(z) = <x 2 (l — a)(l — 2a) — 2(1 — a)kz + 2( a 1 n T 0 1 ) then it is easy to 

2 ‘2 — ‘2 cl 

see that for any M > 0 there exists a cm such that b(z) > M — c m z. We construct 
now the following CIR process 

ft=x o + [ (M - c M f s )ds + 2(1 - a)cr / \JTsdW s . 

Jo Jo 

Using a comparison theorem for stochastic differential equations (see [E], prop. 5.2.18) 
we deduce that z t > ft > 0 (choosing big enough M > 0 in order f t to be strictly posi¬ 
tive and also the inverse moment bound holds for f t ) and using the inverse exponential 
moments of [9] we have the desired result. □ 


Theorem 1 Under Assumption A we have 

E|ay — y t \ 2 < CA 2a{a ~ 1/2 \ 


and therefore the order of convergence is at least a(a — 1/2). 

Proof. Using Ito’s formula on \x p — y p \ 2 for some stopping time p, we obtain 


E\x p -h p f = / 2kE(x s - h s )(x s - y s ) + a 2 E(a:“ - y a s sgn{z s )) 2 ds 


< 


(2kE\x s — h s | 2 + 2kE\x s — h s \\h s — y$\ + a 2 E(x“ — ygSgn(z s )) 2 ds. 


But 


E« - y a s sgn(z s )) 2 < 2E(x“ - y a s ) 2 + 2Ey 2a (l - sgn{z s )) 2 

= 2E(^ - y a s ) 2 + 8P(^ < 0)E (y 2a \(z t < 0)) 

A a-i/2 

< 2E(x“ — y a s ) 2 + C - d — 

g £^2a — 1 

< CA 2a{a ~ l/2) + 2E(x“ - h a s ) 2 . 
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We have used the fact that 

E\y a s - K \ 2 < CE\y s , - h s \ 2a < C (E\y s - h s | 2 )“ < CA 2a{a ~ l/2) . 
Moreover, by Young inequality, 

E\x s — h s \\h s — ys\ < l/2E|x s — h s \ 2 + l/2E|/z s — y- s \ 2 < CA 2a ~ l + l/2E|x s — h 
Therefore, 

E|x p - h p | 2 < CA 2 ^- 1 ^ + f P (3fcE|z s - h ,| 2 + 2a 2 E(a:“ - h a s ) 2 ) ds 

Jo 

Setting now ds and using the inequality 


we have 


W +v. 

x a — y a \\x 1 ~ a + y l ~ a | < 2a\x — y\ 


E\x p — h p \ 2 < CA 2a ( a + / E|x s — y s \ 2 (3ks + y s ) ds 


(4) 


Dehne the stopping time 

t i = inf{s G [0, T] : 3 ks + > l}. 

Now, for p = 77 , we use the change of variables setting u = 3 ks + and therefore 
s = r u obtaining 

E(x Ti - h Tl ) 2 < CA 2 < a ~ 1/2) + [ E\x Tu - h T fdu. 

Jo 

Using Gronwall’s inequality we obtain, 

E\x n - h Tl | 2 < Ce l A 2a{a ~ 1/2) . (5) 

Going back to (4), for p = t G [0, T], we have under the change of variable u = 3 /cs+ 7 s , 


E{x t -h t ) 2 < GA 2a(a - 1/2) + E 


f3fcT+7T 


x Tu ~ K u \ 2 du 


< CA 2a(a_1/2) + / E (l { 7 T > u} |x Tu — h r J 2 ) du. ( 6 ) 

Jo 

Using (5) and noting that 

E ( I {lT > u} \x Tu - h T J 2 ) = P( 7 T > u )E (\x Tu - h Tu \ 2 \{^ T > u)) 


E (x t - h t ) 2 < CA 2a{a ~ l/2) (1 + / P( 7t > u)e u du ) . 


we arrive at 


Note that, 

P(7T > u) < 4e(6 7T ), 
e u 

therefore, using the inverse exponential moments of the true solution we deduce that 

E (x t - h t ) 2 < GA 2a(a - 1/2) 

□ 
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